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Abstract 

We suggest that the atmospheric neutrino anomaly observed in the Super-Kamiokande 
(and other) experiments results from the combined effects of muon-neutrino to tau- 
neutrino oscillations with a Am 2 value of approximately 0.4 eV 2 and oscillations be- 
tween muon neutrinos and electron neutrinos (and vice-versa) with 0.0001< Am 2 < 
0.001 eV 2 . With an appropriate choice of a three-neutrino mixing matrix, such a 
hypothesis is consistent with essentially all neutrino observations. 
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1 Introduction 



Statistically significant anomalies have been reported in observations of so- 
lar neutrinos [1-5], atmospherically-produced neutrinos [6-8], and accelerator- 
generated neutrinos [9]. These observations raise the question whether they 
can all be explained by a simple three- neutrino oscillation model [10-15]. So 
far, data on the atmospheric neutrino anomaly have been fitted only to a two- 
neutrino oscillation hypothesis. We believe that this may lead to an incorrect 
conclusion regarding the oscillation process and its associated parameters. 

Several aspects of the Super-Kamiokande data motivate us to propose that 
muon-neutrino to tau-neutrino oscillations are observed with Am 2 around 
0.4 eV 2 and that, furthermore, oscillations between muon neutrinos and elec- 
tron neutrinos (and vice-versa) are seen with 0.0001 < Am 2 < 0.001 eV 2 . 
First, low-energy (< E >=0.75 GeV) "overhead" (L~20 km) events show a 
highly significant deficit of muons, relative to electrons, indicating an oscilla- 
tion process with Am 2 >0.04 eV 2 . This would be consistent with the positive 
oscillation result obtained in the LSND experiment which, when combined 
with negative results from other experiments, requires Am 2 ~ 0.4 eV 2 . Sec- 
ond, Super-Kamiokande observes an additional, significant muon deficit for 
"upward" events for which the pathlength L~10,000 km. These events also 
show some enhancement of electron events, suggestive of muon-neutrino to 
electron-neutrino oscillations. 



2 Assumed Model for Neutrino Oscillations 



In what follows, we assume the simplest general case of three-neutrino mixing: 
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where v e ,v^,v T are, respectively, the electron, muon, and tau neutrino states, 
and i>i,i>2, and 1*3 are the neutrino mass eigenstates with masses mi < m 2 < 
m 3 . For simplicity, the unitary mixing matrix U is assumed to be real and 
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expressed as: 
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where c stands for cosine and s for sine of the appropriate three angles, 61,62, 
and 6 3 . 

Oscillations of neutrinos of type a into neutrinos of type (5 occur with a 
probability given by: 
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where = 1 if a = (3 and 5 a p = if a 7^ L is the length of the neutrino 
flight path, and E is the neutrino energy. If m is given in eV, E in GeV, and 
L in km, then the factor of L/4E in the sine term becomes 1.27 L/E. 

We now discuss the interpretation of various observations in terms of this 
simple model. Given our hypothesis regarding the atmospheric neutrino data, 
we shall assume that m 3 3> m 2 > mi. We define: 
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With M 2 > 777 2 . 

For the mixing parameters Uij obtained in our analysis, the value of M 2 
allowed by LSND is highly constrained by the results [16,17] of the Bugey 
reactor experiment (M 2 > 0.2 eV 2 ) and of the CDHS oscillation search 
(M 2 < 0.4 eV 2 ). We shall assume that M 2 = 0.4 eV 2 for reasons noted 
below. To produce a significant effect in Super-Kamiokande the value of m 2 
must be at the 0.0001 eV 2 level or larger. The CHOOZ reactor experiment 
[19] has shown that Am 2 is smaller than 0.001 eV 2 . We shall therefore require 
that 0.0001< 777 2 < 0.001 eV 2 . 
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3 LSND 



The LSND probability to observe oscillations of muon neutrinos into electron 
neutrinos is given by: 



P(LSND) = -AU 23 U 13 U 21 U ±1 sin 2 
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In the LSND experiment [9] L=30 m and 36 < E < 60 MeV, so that for 
m 2 < 0.001 eV 2 the third term in the above expression can be neglected. The 
first two terms can be combined by approximating M 2 + m 2 R3 M 2 and using 
the unitarity of the mixing matrix: 

P(LSND) = A[U 23 U 13 ] 2 sin 2 (^) = A(s 2 c 3 s 3 ) 2 sin 2 (^) (6) 



The measured value is P(LSND)=0.0031 ±0.0011(stat) ±0.0005(syst). When 
other neutrino data are also considered, the smallness of P(LSND) arises in 
part from the smallness of S3. For M 2 = 0.4 eV 2 , L=30 m, and a mean value 
of E=42 MeV, we obtain s 2 c 3 s 3 =0.0784. 



4 Atmospheric Neutrinos 



Super-Kamiokande will observe oscillations of muon neutrinos into tau neu- 
trinos with a probability given by: 
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For the oscillation parameters considered in this paper, the third term can be 
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neglected and the first two terms combined to obtain: 

M 2 L 

P(SK; V(l -> v T ) = 4[s 2 c 2 c 2 ] 2 sin 2 — — = 2[s 2 c 2 c 2 ] 2 (8) 



where the variation in E and L allow us to average the sine-squared factor to 
a value of 0.5. 

The smallness of the LSND effect and the failure to observe electron-neutrino 
disappearance in reactor experiments [16,19] indicate that the overall deficit of 
muons observed in the Super-Kamiokande atmospheric data for zenith angles 
satisfying cos9 z > —0.6 results primarily from oscillations of muon neutrinos 
into tau neutrinos. Using the result that P(SK)=0.30 for such events [6], we 
obtain s 2 c 2 c§ = 0.387. Combining this with the LSND results that S2C3S3 = 
0.0784, we find 9 2 = 26.5° and 9 3 = 10.3°: 

We quote results using central values for various parameters. Clearly, the un- 
certainty in an input such as P(LSND) will produce uncertainties in a param- 
eter such 9 3 (of order ±2°). Our purpose is to show that these central values 
can account for essentially all observations. 

Super-Kamiokande will also observe oscillations between muon neutrinos and 
electron neutrinos with a probability given by: 
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(M 2 + m 2 )L] 



-4U23U13U22U12 sin" 



-4U22U12U21UU sin" 



AE 
M 2 L 

AE 
m 2 L 
~A~E 



(9) 



The sine-squared factors of the first two terms of this expression each average 
to a value of 0.5 and the two terms can be combined to yield: 

P(SK;v^ <- v e ) = 2[f/ 23 f/ 13 ] 2 - 4C/ 22 C/ 12 C/ 21 C/nsin 2 (^) (10) 



For our values of 9 2 and 9 3 , the first term is small (~0.02), justifying our 
assumption that, for most of the zenith-angle range, the atmospheric muon 
anomaly is caused by muon-neutrino to tau-neutrino oscillations. However, 
the second term will produce observable effects for "upward" going events for 
which is assumed to be of order one or greater. We return to this point 
in the discussion below. 
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We have examined the question of matter (MSW) effects on oscillations be- 
tween atmospheric electrons and muon neutrinos [18]. For the large mixing 
angles obtained for our model, resonance effects do not have a significant im- 
pact on our analysis. However, for m 2 values near 0.0001 eV 2 , oscillations 
between electron and muon neutrinos are strongly damped for neutrino en- 
ergies above 1 GeV. Such damping effects decrease rapidly with increasing 
values of m 2 , and are insignificant at the upper end (0.001 eV 2 ) of the m 2 
range considered here. 



5 Solar Neutrinos 



The MSW mechanism is not expected to have a significant impact on solar 
neutrinos for the large mass-squared differences considered in this paper, so 
that the probability for a solar electron neutrino to remain as such is: 



P(solar; v e — > v e ) = 1 - A[U n Un} 2 sin 2 
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For the E and L of solar neutrinos detected on Earth, the sine-squared terms 
in the above expression average each to 0.5 and: 

P(solar; v e -> v e ) = 1 - 2(s 3 c 3 ) 2 - 2(s lCl c 2 ) 2 (12) 



Within the uncertainties of solar neutrino theory, essentially all solar neutrino 
observations (with the possible exception of those from the Homestake experi- 
ment) are consistent with P(solar)=0.5 [20]. Using this value and our previous 
result for angles 9 2 and 9 3 , we obtain 9i = 37.6°. 



6 Reactor Experiments 



The general expression derived for solar neutrinos also applies to reactor ex- 
periments. Such experiments have not been sensitive to m 2 values less than 
0.001 eV 2 . However, oscillations with M 2 = 0.4 eV 2 are, in principle, observ- 
able. For values of E and L that allow an averaging of the oscillation factor, 
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we expect: 

P(reactor; v e — >• v e ) — 1 — 2(s 3 c 3 ) 2 



(13) 



With # 3 = 10.3°, we predict P(reactor) = 0.94 compared to measured values 
of 0.99±0.01(stat)±0.05(syst) for Bugey [16] and 0.98 ±0.04(stat) ±0.04(syst) 
for CHOOZ [19]. Given the size of the errors, this is reasonable if not perfect 
agreement. In the Bugey experiment, the detectors were sufficiently close to the 
reactor to be sensitive, in principle, to modulations in the positron spectrum 
produced by the sin 2 ^rr- factor. We discuss this below. 



7 Summary and Discussion 



We have found that most features of oscillation-related neutrino data can 
be explained by the following mixing matrix, corresponding to the angles 
0i = 37.6°, 9 2 = 26.5°, and 9 3 = 10.3°. 



U 
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-0.61 

0.15 



0.60 0.18 
0.66 0.44 
- 0.45 0.88 



(14) 



and by assuming that neutrino masses m 1 ,m 2 ,m 3 satisfy m 3 >> m 2 > mi, 
that 0.0001 < (m 2 - mj) < 0.001 eV 2 , and that M 2 = m 2 - m 2 = 0.4 eV 2 . 
This choice of parameters appears to be in reasonable accord with essentially 
all neutrino observations, including the zenith-angle behavior of the Super- 
Kamiokande atmospheric neutrino data for "upward" going (cos# 2 < —0.6) 
events. We demonstrate this by assuming that L is sufficiently large that all 
sine-squared factors involving mass-squared differences average to 0.5. We find 
that: 



= P{SK- *v - Vp) = f4 + f4 + f/ 2 4 3 
P ee = P(SK; v e - v e ) = Uti + Ut 2 + Uf 3 

P e , = P,e = P(SK;v^ <-> v e ) = 2[U 23 U 13 ] 2 - 2U 22 Ui 2 U 2 iUu (15) 



The double ratio R for these "upward" going events: 

(Nft/Ne) measured 



(Nn/N e ) no oscillation 1 ; 
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can be estimated to be: 



R = (2P^ + P e ,)/(P ee + 2P, e ) 



given the fact that atmospheric are produced at approximately twice the 
rate of v e . We obtain the result that R=0.44 compared to the measured value 
of 0.41±0.04 at the largest zenith angles. 

We show in Fig. 1 the predicted ratio R as a function of L/E for three values 
of m 2 = m| — ml- The curve of R vs. L/E has been smeared by a resolution 
function that is Gaussian in ln(L/E) with a width that corresponds to a factor 
of three uncertainty in the L/E measured for any given event. In an exper- 
iment such as Super-Kamiokande the angle between the incoming neutrino 
and the direction of the detected lepton can be quite large [6], leading to a 
significant uncertainty in L. The general shape of the curve of R vs. L/E in 
Fig. 1 depends only weakly on the assumed smearing, with factors of two or 
four giving qualitatively the same results as our assumed factor of three. We 
have also ploted in Fig. 1 the latest Super-Kamiokande results [6] for which a 
reasonably good fit is obtained for m 2 values around 0.0003 eV 2 . 

The most recent Super-Kamiokande paper [6] makes available the ratio of the 
number of measured to the number of expected events, assuming no oscil- 
lations, for electrons and muons separately. We show these data in Fig. 2, 
with the overall flux normalization reduced by 12% to give a best fit to the 
three-neutrino mixing model presented in this paper. Such a normalization 
shift is well within the errors of the theoretical flux assumed by the Super- 
Kamiokande Collaboration. Our model does not give a very good fit to the 
individual electron and muon ratios. However, it is important to note that 
the double ratio R shown in Fig. 1, for which we obtain a good fit, is much 
less sensitive to systematic errors than are the individual electron and muon 
ratios. Until these systematic errors are reduced by improved measurements 
of atmospheric cosmic-ray particle production, one should keep an open mind 
about the significance of individual distributions in the atmospheric neutrino 
data. 

As with previous attempts to explain all neutrino anomalies with a simple 
three-neutrino mixing model, the one outlined in this paper is subject to the 
resolution of various experimental inconsistencies and ambiguities. For ex- 
ample, we have assumed that the LSND result is an oscillation effect rather 
than unexplained background. It must be noted that the KARMEN experi- 
ment has not confirmed the LSND result. For our assumed parameters, recent 
KARMEN data [21] should have yielded approximately two signal events in 
addition to roughly three background events, whereas no oscillation candidate 
events were seen. However, the LSND collaboration has reported additional 
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Fig. 1. The predicted double ratio R as a function of L/E for the mixing matrix 
with 6>i = 37.6°, 2 = 26.5°, 3 = 10.3° and for three values of m 2 = m\- m\. The 
value of M 2 = m 2 — to 2 , = 0.4 eV 2 . The data are the most recent Super-Kamiokande 
results [6]. 

evidence for oscillations in events in which the neutrino was produced by de- 
cays in flight [22]. 

With respect to the atmospheric neutrino data of Super-Kamiokande, our 



model with Q x = 37.6°, 2 = 26.5°, 3 = 10.3°, M 2 



m% 



mr, 



0.4 eV 2 



and m 2 = ml — ml = 0.0003 eV 2 gives a good overall fit to the data. To 
distinguish our model from one with a single oscillation process with a small 
Am 2 (~ 0.0022 eV 2 ) requires precise measurements of multi-GeV, "overhead" 
(cos0 2 ~ 1) events. Our model predicts R~0.72 for such events whereas the 
small- Am 2 model yields R~l. Present Super-Kamiokande data are inconclu- 
sive in this low-L/E region. 
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Fig. 2. Separate Super-Kamiokande electron and muon data with overall flux nor- 
malization reduced by 12% for a best fit to the three-neutrino mixing model pre- 
sented in this paper with m 2 = 0.0003 eV 2 . 

Our model requires that solar neutrino deficits observed on Earth should show 
no energy dependence and no MSW effects. The data from the gallium experi- 
ments (GALLEX and SAGE) and from Super-Kamiokande, which sample very 
different parts of the solar neutrino energy spectrum, are in general agreement 
with this prediction. However, the Homestake Cl-37 experiment has histori- 
cally reported a substantially larger solar neutrino deficit [23] than is observed 
in these other experiments. We note, however, that the measured Homestake 
solar neutrino flux of 2.56 ±0.16(stat) ±0.16(syst) SNU represents 0.403 ± 
0.025 ±0.025 of the expected flux computed by Turck-Chieze and Lopes [24], 
not far from the value 0.50 assumed in this paper. A final interpretation of 
all the data awaits a better understanding of both the theory and of possible 
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systematic errors in various experimental results. 

Our model agrees within one standard deviation with the fluxes measured in 
reactor experiments. However, for 9 5 = 10.3°, the Bugey results [16] would 
be expected to show some modulation (at the level of several percent) of the 
shape of the measured positron spectrum. We note, however, that statistical 
and systematic uncertainties in the LSND result permit values of 83 as small 
as 8°, which would reduce such modulations below the level of detectability. 
Our choice of M 2 = 0.4 eV 2 , the maximum allowed by the CDHS oscillation 
search [17], is dictated by the requirement that 63 be small. 

We note that our model predicts large effects for terrestrial, long-baseline 
experiments such as K2K, MINOS, and those planned for Gran Sasso. Fur- 
thermore, if muon neutrinos oscillate into tau neutrinos with M 2 = 0.4 eV 2 , 
then observable effects are just beyond the reach of the present short-baseline 
experiments CHORUS and NOMAD but accessible to a new, improved exper- 
iment such as proposed by the TOSCA collaboration. 

Finally, we remark that our model differs from previous versions of simple 
three-neutrino mixing [10-15] in its unique combination of a M 2 value that 
can explain the LSND results and a m 2 value significantly larger than is usually 
assumed in solar neutrino models. It is, of course, possible that eventually no 
simple three- neutrino mixing model will give a satisfactory account of past, 
present and future experimental results. In that case, more complicated models 
will have to be considered [25]. 
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